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Abstract 

We review current understanding of the non-equilibrium dynamics of collective quan- 
tum systems. We describe an approach based on the Hamiltonian formulation of su- 
perfluid hydrodynamics. It is shown that, in the presence of constant energy pumping, 
the nonlinear coupling of fluctuations in the density and entropy strongly affects the 
noncquilibrium dynamics of the system. We use the results obtained to analyze the 
properties of out-of-equilibrium superfluid 4 He and of exciton polariton Bose-Einstein 
condensates, both in semiconductor quantum wells and in graphene layers in presence 
of high magnetic field. 



1 Introduction 

Superfluidity is one of the most fascinating macroscopic quantum effects in condensed mat- 
ter. It exhibits itself as a frictionless, collective movement of fluids - either liquids or gases. 
After being discovered in 1938 by Kapitza [1] and by Allen and Misener [2] in their famous 
experiments on helium 4 He (the stable isotope of helium, of atomic mass equal to 4) at 
temperatures below ~ 2 K, this phenomenon exerted an influence on the development of 
condensed matter physics for the rest of the 20 th century. It is worthy of note that su- 
perconductivity (that is, a complete loss of electrical resistivity at low temperatures) was 
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later explained as the superfluid motion of electronic Cooper pairs in metals [3,4]. We also 
mention that the apparent formation of a superfluid fraction in solids, called supersolid, 
has been recently observed at mK temperatures [5]. The fundamentals of the theory of 
superfluid liquids were developed in the pioneering works of Lev Landau [3] and Isaac Kha- 
latnikov [6]. It was subsequently established that superfluidity can be understood in terms 
of Bose-Einstein condensation (BEC) of particles and establishment of a long-range order 
in the system, known as quantum coherence [7, 8]. In particular, atoms in the ground state 
of a supefluid are described by a single macroscopic wave function ^ = y/pexp(i6), where 
p is the condensate density and 9 is the phase. 

In the last few years, interest in superfluid systems has further increased because of 
the discovery of Bose-Einstein condensation in atomic gases at pK temperatures [9, 10]. 
For the discovery and study of this phenomenon Cornell, Ketterle and Wieman shared the 
2001 Nobel Prize in Physics. Furthermore, it was recently found that excitons, which are 
bound electron-hole pairs in semiconductor nanostructures, can also form a state where 
macroscopic quantum coherence strongly affects their dynamics [11-13]. Other examples 
of systems in which quantum coherence is important are provided by polaritons, which 
are a resonant coupling of excitons with cavity photons in semiconductors [14-16], and by 
magnitopolaritons, which are generated in graphene embedded in a microcavity when mag- 
nitoexcitons are coupled with photons in a strong magnetic field [17]. Geim and Novoselov 
have received the 2010 Nobel Prize in Physics for the discovery and study of graphene. 

It was found that, under certain conditions, the polaritons form a superfluid liquid that 
is in many aspects similar to superfluid 4 He [16]. Because of the extremely small polariton 
effective mass, the latter can exist as a superfluid at temperatures much higher than that 
of the temperature of normal-to-superfluid transition in liquid 4 He. The example of BEC 
and superfluidity of magnitoexcitons in graphene in strong magnetic fields is of special 
importance because graphene is considered as a hot candidate for future applications in 
electronics. 

The key feature of the exciton polariton superfluid in microcavities is that it can only 
form under nonequilibrium conditions: the composite particles have a finite lifetime and the 
particle losses must be compensated by continuous production (typically, resonance pump- 
ing by laser radiation is used [18]). Because of that, there is currently a strong interest in ki- 
netic phenomena in non-equilibrium superfluid liquids [19-21]. Although out-of-equilibrium 
superfluids can be studied, in principle, in a number of systems, superfluid helium 4 He in 
a macroscopic cavity presents a useful testbed for experimental studies of general kinetic 
and transport properties. This potentially fruitful approach has already yielded important 
results [22-24]. Specifically, it has been demonstrated that the main mechanisms of energy 
relaxation in the superfluid liquid is governed by the nonlinear dynamics of second sound 
oscillations and its coupling with ordinary sound, i.e. with density waves in the bulk super- 
fluid. (Second sound in a quantum system is a slightly dissipative entropy wave that can be 
understood as a density wave in the gas of quasiparticles [6]: see Sec. 2). Non-equilibrium 
conditions were achieved in these experiments through continuous driving at a resonant 
frequency of the cavity by a periodic drive applied to a heater immersed in the superfluid 
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[22]. It was found that, because of strong nonlinearity of the second sound interactions, 
the perturbations play an important role in the non-equilibrium dynamics of the system; 
this interaction results in the development of energy cascades and the formation of constant 
fluxes of conserved quantities (e.g. energy and wave action) through the frequency scales 
(cf. Sec. 2.2). 

In this paper, we review recent achievements in the theory of non-equilibrium superfluid 
systems. In Sec. 2, we develop a general formalism that describes its nonlinear dynamics. 
Through this technique, we find steady-state non-equilibrium solutions of the proper kinetic 
equations describing the correlation functions for the density and entropy fluctuations. 
Here, we take advantage of the Hamiltomian approach to describe of nonlinear dynamics 
of non-equilibrium superfluid liquids. In particular, we utilize an approach based on the 
kinetic equation for the correlation functions for the density and entropy fluctuations which 
is a generalization of that used in consideration of classical nonlinear turbulent systems. 
Through analytical computations and numerical simulations, we describe the concept of 
second sound turbulence, which was observed in recent experiments [22]. In Sec. 3, we 
focus on application of the Hamiltonian formalism to a particular system of exciton polariton 
supefluids in microcavities. In doing so, we demonstrate the transition to a superfluid state 
in this system and describe the main characteristics of the Bose-Einstein condensate fraction 
as functions of the interaction parameters of exciton polaritons with each other and with 
an external trapping field. 



2 Nonlinear phenomena in superfluid liquid 

2.1 Hamiltonian formulation of superfluid hydrodynamics in quantum 
liquids 

2.1.1 Superfluid hydrodynamics 

As mentioned above in the Introduction, the Hamiltonian formalism plays a central role 
in the subsequent consideration of kinetic processes in quantum superfluid systems. For 
the reader's convenience, we now sketch the Hamiltonian formulation of the superfluid 
hydrodynamics equations [22,25-29]. The Landau's well-known hydrodynamics equations 
for a quantum liquid that has the total density p and the densities p s and p n for its superfluid 
and normal components are [6] 



dp n . . 9S .„ . 

-g + div j = 0, — + div (Sv n ) = 

dv s „ ( 1 ,\ dji 9Hai. 



(1) 



Here S is the entropy per unit volume, v s and v n are the velocities of superfluid and normal 
motion, j = p s v s + p n v n is the mass flux, U ik = p5ik + p s v s iv s k + Pn v n i v n k is the tensor 
of the momentum flux, \i = de/dp is the chemical potential of the unit mass, and e is the 
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energy per unit volume. Equations (1) are the continuity equations for the density, entropy 
and momentum of the fluid, and the equation for the superfluid velocity, which is the phase 
gradient of the condensate wave function, v s = KVO/m^, where fi is Planck's constant and 
777,4 is the mass of a 4 He atom. The damping terms in Eqs. (1) are omitted. In this review, 
we restrict our consideration to processes for which the thermal expansion of the fluid is 
not an essential feature, and the terms proportional to the thermal expansion coefficient 
kt = (T/p)(dp/dT) of the fluid are therefore neglected. For superfluid 4 He, this coefficient 
is indeed small, with kt ~ 10~ 2 over the range of temperature T < 2.15 K. 

In a non-stationary state, the density, entropy and velocities of the superfluid are not 
equal to their equilibrium values. For an infinitely small deviation from equilibrium, Eqs. 
(1) are reduced to two linear wave equations for perturbations of the density 5p, and the 
entropy per unit mass 5a [6], 



a = S/p, and c is the heat capacity per unit mass. The linear equations (2) describe the 
propagation of two types of waves in the bulk superfluid: density waves in which 5a = 
at velocity u\, and entropy waves in which 5p = at velocity u 2 . These waves are known 
as first sound and second soimd,respectively [6]. In second sound, the fluid density remains 
constant but the temperature varies with coordinates and time as 5T = (dT / da)5a [30]. 
The possibility of propagating of temperature waves that are only weakly dissipative is a 
hallmark of quantum systems, and is in sharp contrast to diffusive heat transfer mechanisms 
in classical media where the temperature waves is damped over a distance of the order of 
its own wavelength [31]. 

2.1.2 Hamiltonian representation of hydrodynamic equations for the super- 
fluid liquid in a cavity 

As we have seen above, the dynamics of a superfluid liquid can be described, to a first 
approximation, as the propagation of the density and entropy waves. In what follows, we 
disregard modes that are localized at the fluid boundaries. (For a detailed consideration 
of nonlinear phenomena at the surface of quantum liquids, see Refs. [26, 29, 32]). At finite 
wave amplitude, however, nonlinear effects in the wave interaction play an important role. 
In particular, second sound is characterized by rather strong nonlinear properties [33,34]. 
This can lead to the formation of shock waves (temperature discontinuities) during the 
propagation of finite amplitude waves at short distances from the source [35-39]. The 
velocity of a traveling second sound wave depends on its amplitude and, to a leading order 
approximation, can be written as u 2 = u 2 q(1 + ST), where 5T is the wave amplitude, u 2 o is 
the velocity of a wave of an infinitely small amplitude, and 




d 2 5a 

= ujASp, = u 2 2 A5a. (2) 



Here u\ = (dp/ 'dp) 1 / 2 and u 2 = (p s a 2 T / Pnc) 1 / 2 are the first and second sound velocities, 




(3) 
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is the nonlinearity coefficient. In He-II, the superfiuid phase of liquid 4 He, the nonlinearity 
coefficient 02 of second sound may be either positive and negative, depending on the tem- 
perature and pressure [6,33,39,40]; under the saturated vapor pressure, the nonlinearity 
coefficient is positive («2 > 0) at temperatures T < T a = 1.88 K (like the nonlinearity of 
conventional sound waves in ordinary media), but it is negative in the range T a < T < T\. 
Here, T\ = 2.176 K is the temperature of the superfluid-to-normal transition in bulk 4 He. 

To describe nonlinear effects, we now introduce the Hamiltonian variables for the system 
(1). An arbitrary flow of superfiuid can be described by three pairs of conjugate variables 
(a,p), S) and (7, /) [25, 26]. Here a = h9/m is the superfiuid velocity potential, f3 is the 
phase variable conjugated to S, and 7 and / are Clebsch variables. The Clebsch variables 
define the vorticity in the system; we focus below on the potential motion of the superfiuid, 
so we set 7 = / = 0. The Hamiltonian function of the superfiuid system is given by the 
total energy expressed in the stationary frame of reference, 



H = d 6 r 



Pn 2 1 Ps 2 1 / a \ 

y^ + y« s +£(p, S,p) 



(4) 



Here e(p,S,p) is the energy per unit volume of the superfiuid component in the reference 
frame moving with the velocity v s , and p is the momentum of relative motion of the normal 
component. The total mass flux is expressed through these variables as j = pv s + p = 
pVa + The equations of superfiuid hydrodynamics (1) for the Hamiltonian function 

(4) are derived in terms of the conjugate Hamiltonian variables as follows [25, 26] 

. SH 8H . 5H • 5H 

In Eqs. (5) 5 denotes the variational derivative [43]. In this representation, the momentum 
of the relative motion of the normal component is equal to p = SV/3. 

If the sound waves propagate in an unrestricted superfiuid system, the conjugate vari- 
ables can be expressed in terms of normal coordinates - the normalized amplitudes of the 
running first and second sound waves [25]. In a cavity, the normal modes satisfy appropriate 
boundary conditions. In the important case of a high-quality resonator [22,23,36,44,45], 
the corresponding normal modes are standing waves of first and second sound whose fre- 
quencies are equal to the resonant frequencies, fifc = ui|fc n |, lo^ = U2\k n \, where k n is the 
nth resonant wave vector. 

Variations of the Hamiltonian variables in a wave with respect to their equilibrium values 
are expressed via the sound amplitudes as [25, 26, 41] 

t) = Yl a ¥k(r)(a k - a* k ), 5p(r, t) = ^ p^>k{r){a k + a%), 

k k ( 6 ) 

P(r,t) =Y J hk{r){b k - b%), 6S(r,t) = ^5^ fc (r)(6 fc + b* k ), 

k k 
where the normal coordinates a k and b k are the complex amplitudes of the first and second 
sound waves, respectively, p, S, a and (3 are the normalizing factors, and a star stands 
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for the complex conjugate. The complex amplitudes are utilized to take the wave phases 
into account. The spatial basic functions depend on the geometry of the cavity and the 
boundary conditions. For the simplest case of a high quality rectangular resonator, they 
are (pk(f) = cos(k x x)cos(k y y)cos(k z z). The equations of motion (1) are presented by the 
Hamiltonian equations for the sound amplitudes as follows: 

** = f| - ^a k + if \ ib k = fg - rrS\ + if\ (7) 

(i) (i) 

The dissipation 7^ , and the interaction with an external driving force F k , are included 
phenomenologically [25, 26]; the superscript % = 1, 2 labels the first and second sound modes, 
respectively. Note that this representation provides a quasiclassical limit for the equations 
of motion of a superfluid system that can be used provided that the occupation numbers of 
the corresponding states are sufficiently large, i.e. |afc| 2 , |^fc| 2 3> ^- However, for the purely 
quantum case \dk\ 2 , \bk\ 2 ~ fr, the normal variables must be considered as operators in a 
proper Hilbert space; we will consider examples of such systems in Sec. 3. 

The nonlinear dynamics of the system (7) can be considered in terms of an expansion of 
the Hamiltonian H in a Taylor series over the wave amplitudes a k and H = H2+H3 + . . .. 
The term ii2 is a quadratic function of the amplitudes, H2 = J2k {^k\ a k\ 2 + ^k\bk\ 2 ) , 
that describes the propagation of linear waves [46]. Nonlinear effects such as parametric 
generation and decay of waves correspond to anharmonic terms in the Hamiltonian, with 
the coefficients in the expansion over a& and being the amplitudes of the processes. 
Three- wave processes include (a) the decay of the first sound into two second sound waves, 
(b) Cerenkov emission of the second sound wave by the first sound wave, and (c) inner 
decay and confluence of waves that belong to the same wave mode. As we show in Sees. 
2.2 and 2.3, these processes are responsible for the generation of the direct and inverse 
energy cascades in frequency space for a superfluid liquid. In actual experimental regimes, 
one of these processes may sometimes dominate (see Refs. [22, 23, 49-52] and discussions 
in Sec. 2.3 below). In what follows, we disregard interactions which involve four or more 
waves. Therefore, the higher order terms in the expansion of the Hamiltonian over the wave 
amplitudes are omitted. These higher-order processes are responsible for the isotropisation 
in the system [47, 48] and should be taken into account in a more general theory. 



2.1.3 Statistical description 

Under experimental conditions [22, 23, 41, 44, 53-56], many degrees of freedom (that is, wave 
modes) in a superfluid liquid are excited simultaneously and interact with each other. In 
this far-from-equilibrium state of a nonlinear system, the dynamics should be described 
statistically. The theory of wave turbulence, initially formulated for classical wave fields 
[46, 57], presents a useful framework for this description. Below, we extend this approach 
to the case of nonlinear quantum systems and, in particular, we study kinetic phenomena 
in a superfluid system described by the three- wave Hamiltonian equations (4)— (7). In this 
statistical approach, the state of the quantum system is characterized by simultaneous 
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correlation functions of the normal variables, averaged over the ensemble of waves 



(al(t)a kl (t)} = N k (t)5(k - k x ), {b* h (t)b kl (t)) = n k (t)S(k - k x ). (8) 

The correlation functions N k (t) and n k (t) are the "occupation numbers" of the first and 
second sound wave modes. The time evolution of the correlation functions obeys kinetic 
equations that are similar to the Boltzmann equation [59], 

^ = h[N k ,n k ], ^ = h[N k ,n k ], (9) 
with h[N k ,n k ] and l2[N k ,n k ] being the "collision" integrals [42,46,49,60]. 

2.2 Kinetic and nonstationary phenomena in superfluids 
2.2.1 Equilibrium fluctuations 

The kinetic Eqs. (9) can be used for studying the energy relaxation processes in a superfluid. 
In Sec. 2.2, we focus on the case where one of the parametric processes - the decay of first 
sound to two second sound waves and the inverse process of confluence of two second sound 
waves into a first sound wave - dominates. For superfluid 4 He there is a large window 
of bath temperatures and wave amplitudes, in which this process is relevant [25,51]. The 
opposite situation where the mutual transformation of the first and second sound waves is 
negligible, and the wave dynamics is governed by the inner nonlinearity of the wave modes, 
is described in Sec. 2.3. 

In the absence of an external energy input, the system tends towards thermodynamic 
equilibrium where only temperature fluctuations exist. Such an equilibrium state is de- 
scribed by a general Rayleigh-Jeans distribution [58] 

T T 
iifc + (vo«) + \l 0J k + (Vofe) + jU 

which are the stationary solutions of the kinetic equations (9). The vector vo is the constant 
drift velocity of the excitations in the superfluid with respect to the stationary frame of 
reference, and T is the temperature. The term (vgfc) is simply the Galilean transformation 
to a moving frame of reference. In Eqs. (10) /i, fi' play a role of the effective chemical 
potentials for the quasiparticles, which obey the relation \x = 2//. Let us note that (i and 
fj! do not coincide with the chemical potential of atoms in the medium; the latter always 
tends to zero for a Bose-Einstein condensed macroscopic superfluid system. For the fluid 
at rest, one has vo = 0. In this case, in the low-frequency limit, Q< k , u k <C /i, //, one finds 
from Eqs. (10) that N k = n k /2 = T/fi. Thus, in equilibrium, the long-wave distributions 
do not depend on k and are proportional to the temperature T of the system. In the high 
frequency limit, U k , uj k S> /i, //, one obtains N k , n k oc k~ l . 
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2.2.2 Direct and inverse cascades 



Next we turn our attention to the energy transfer mechanisms in a non-equilibrium super- 
fluid in presence of a continuous energy input. In experiments with superfluid helium, this 
corresponds to the case where the system is continuously driven by an external periodic 
source of heat or pressure [22, 23, 51, 52]. The energy is absorbed in the superfluid through 
bulk dissipation and viscous drag at the resonator walls. For the polariton Bose-Einstein 
condensates the role of the energy pump is played by laser radiation. In this case, the 
effective damping is given by a finite lifetime of the photons in an optical cavity, see Sec. 3 
for a discussion. In both examples, the external driving force is periodic, that is, localized 
in frequency space (we denote the latter as arrive)- 

The steady, non-equilibrium state of the system is described by the equations 

h[N k ,n k }=0, I 2 [N k ,n k }=0. (11) 

An important property of Eq. (11) is that the integrals I\ and I2 contain only power-law 
functions [42]. Because of that, the steady-state solutions of Eqs. (11) acquire a power-law- 
like form 

N k = Ak s , n k = Bk s . (12) 

The scaling indices in Eqs. (12) are equal to s = —9/2 and s = —4 [42]. These solutions 
describe the propagation of constant fluxes of two different conserved quantities, or the 
"integrals of motion" of Eqs. (9), through the scales. It can be shown that these are the 
wave energy, £ = J dk [Vl k N k + uj k n k ], and the wave action Af = J dk [2N k + n k ] (the 
last is sometimes also referred to as "the number of waves" [46]). To show that, we can 
rewrite the kinetic equations (9) in the form of the continuity equations [42] 

d£ dP n dAf dQ . , 

^ + ^ = °< -dT + ^ = °- (13) 

From here it follows that both £ and Af are conserved over time. 

The fluxes P and Q describe the propagation of the corresponding integrals of motion in 
K-space. In the steady state, where d£/dt = dAf /dt = 0, the fluxes P and Q are constant 
in accordance with Eqs. (13). For the non-equilibrium distributions (12), the energy flux 
is equal to P = v\C\p~ l ABy^^ 2 at s = —9/2, and the flux of wave action is equal to 
Q = —V2P~ x ~i\AB at s = —4, where 72 = 1*2/^1 is the second-to-first sound velocity 
ratio, and v\ and V2 are positive dimensionless constants which depend on temperature 
[42]. Thus, in a steady state, the fluxes P and Q have opposite signs. Specifically, the 
energy flux P is positive and therefore, the corresponding integral of motion, i.e. the wave 
energy £, propagates from the driving frequency scale thrive towards the high frequency 
spectral domain, i.e. the direct cascade. As a consequence, the steady-state distribution 
(12) at s = —9/2 forms at frequencies higher than arrive- The Q flux is negative, so the 
wave action, Af, is transferred towards the low-frequency domain and results in formation of 
the inverse cascade. The corresponding solution (12) at s = —4 is established at frequencies 

U) < ^drivc- 
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Note that solutions of the type (12) are similar to the classical turbulent spectra of 
surface waves on deep water [46] or spectra in optical turbulence [61]. In the latter cases, 
however, the inverse cascade appears as a result of the direct four-wave interaction of waves 
that explicitly preserves the number of waves [46]. In a quantum superfluid system, the 
interaction between sound waves is of the three-wave type, and the corresponding integral, 
AT, is conserved as the result of nonlinear interaction between two wave modes - the density 
and entropy waves. Both the direct and inverse cascades have been observed in recent 
experiments with superfluid 4 He in a cavity [23, 24, 55], see Sec. 2.4. 

2.2.3 Propagation of perturbations in A'-space 

We now consider the nonstationary processes through which the cascades are formed. In 
doing so, we use the smallness of the parameter 72 = 112 ju\ to simplify the description of 
the superfluid dynamics. Note that the value of 72 is typically small; for example, in He-II 
one has 72 ~ 0.1 over a wide range of temperatures [6]. 

It follows from the conservation of momentum (fc-vector) and energy (frequency) in 
the decay process that, at small 72, first sound of wave vector k only interacts with the 
second sound whose wave vector is close to k' = k/2^2- This results in locality of the wave 
interaction in A-space. As a consequence, the set of kinetic equations (9) can be reduced to 
a Fokker-Plank differential equation that simplifies the analysis of nonstationary processes 
in the superfluid liquid (cf. Eq. (14)). 

Another peculiarity of the kinetic equations for the decay process is that the charac- 
teristic evolution time for the first sound is very short compared to that for second sound, 
T i/ T 2 ~ 72- Because of this great difference in time scales, we can simplify the whole 
system (9) by considering the first sound waves as a "fast" subsystem, which is always in 
quasi-equilibrium with the "slow" second sound subsystem. Consequently, at large time, 
the wave dynamics is governed by the effective kinetic equation hf. = l2[N^ e \ n^], where 
N [ * e) = n k >/2 is the quasi-equilibrium first sound distribution. Taking into account the 
above-mentioned locality of the interaction, the kinetic integral I2 can be expanded in a 
series over the 72 parameter. To the main approximation, it reads [42] 

1 d 2 An 3 d 2 1 \ 
nk =VW \ k nk WV k )- (14) 

Equation (14) governs the long-time evolution of the system. It is similar to the Fokker- 
Plank kinetic equation for a rarefied gas [59], which describes the diffusion of gas density 
perturbations. In our case, the perturbations of the distributions of coupled first and 
second sound waves diffuse in A-space. However, the r.h.s. of Eq. (14) is the fourth 
order differential operator and hence it corresponds to hyper diffusion. It is worth noting 
that dynamical equations similar to (14) arise in the description of nonlinear light pulse 
propagation in optical waveguides [61] and of long waves on the ocean surface [62, 63]. 
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2.2.4 Formation of the cascades 



Equation (14) provides a useful tool for the analysis of nonstationary processes in super- 
fluid systems. In this Section, we focus on the processes through which the steady-state 
distributions (12) are formed. We suppose that the driving force is applied at t = in a 
step-like manner. Assuming self-similarity of the time evolution of the system, one obtains 
the instantaneous distribution in the high-frequency spectral domain at times t > in the 
following form [45] 



In Eq. (15) we have used the frequency representation with u = 112k to simplify the com- 
parison with experiment, g is a self-similar function, q and p are positive scaling exponents, 
and we denote r = to — t with the finite formation time to ~ T2 where T2 is the second sound 
nonlinear interaction time at oj = arrive- I n accordance with the discussion in Sec. 2.2.3, we 
focus on the second sound distribution n. The first sound distribution can be found from 
an adiabatic relation in each moment of time. The values of the q and p exponents and 
the function g can be found by direct numerical solution of Eq. (14). These results will be 
discussed in more detail elsewhere. 

It follows from Eq. (15) that the formation of the direct cascade in the high-frequency 
spectral domain, u) > Wd r ;vc> is of the "explosion type" [46], with a finite formation time 
~ T2 for the finite capacity spectrum [64]. The transient process can be understood as the 
propagation of a formation front in frequency space from low to high frequencies, whose 
position at time t is ojf(t) oc (to — t)~~ p . 

The formation of the inverse cascade in the low-frequency domain oj < ^drivc can be 
described by the same dependence (15), but with different exponents q and p and with 
r = T2 + t. In experiments with the superfluid liquid in a cavity, the accessible range of 
resonant frequencies is bounded from below by the fundamental frequency of the cavity, 
^0 ~ U2/L. If the driving force is applied at a high resonant frequency, arrive S> the 
build-up process of the inverse cascade requires a time tf ~ ^(^drivc/^o) 1 ^ that is much 
longer than is needed for formation of the direct cascade. This conclusion is in qualitative 
agreement with the experimental observations [23,65]. 

2.3 Numerical simulations 
2.3.1 Second sound turbulence 

To capture the details of the wave cascades' formation processes, we have undertaken a 
numerical study of the dynamics of nonlinear waves in superfluid 4 He within the high-Q 
resonator. In this investigation, we use parameters of the model which correspond to those 
in experiments [22, 23, 45] with superfluid 4 He. The latter were conducted at temperatures 
relatively close to the temperature of the superfluid transition, (T\ — T)/T\ < 0.05. As 
pointed out in Sec. 2.1.2, the nonlinearity coefficient Q2 grows as T — > T\ and hence the 
interaction between the second sound waves becomes important. On the other hand, the 
interaction between second sound and first sound is relatively small in this regime and can 




(15) 
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to a first approximation be neglected. Accordingly, the simulations only take account of 
terms containing the second sound amplitudes bk- 

The frequency of second sound depends on its wave vector k as [66] 

w = u 20 k [l + X fk 2 + ...], (16) 

where £ = £o(l — T/T\)~ 2 ' 3 , £o ~ 2 — 3 A, and Ao ~ 1. The dispersion of second sound 
is significant within a close vicinity of the superfluid transition (i.e. for T\ — T < 1 fiK) 
but is very weak in the temperature range T < 2.1 K. However, the presence of non-zero 
dispersion is of key importance for the wave dynamics in a superfluid liquid. It is well 
established that if the dispersion is exactly equal to zero, duifdk = 0, the main mechanism 
of energy transport through the frequency scales is associated with the creation of a shock 
wave [67]. When the dispersion is positive but small, du2o/dk > (e.g. in the case of 
second sound in He II) , the resonant three- wave interaction is only relevant within a narrow 
cone of k-vectors around the direction of propagation of the wave with, k — k(j r i ve [46,68]. 
The relative phases of waves with k-vectors almost collinear to k<j r i ve (i- e - within the cone) 
are random, and a kinetic equation for waves can be used to describe the cascade-like 
propagation of energy in the wave system [46,69]. However, the interaction of waves with 
non-collinear wave vectors is controlled by higher nonlinear terms and is relatively small. 
That is, the wave distribution in a superfluid is nearly one-dimensional, provided that the 
dispersion given by Eq. (16) is weak. The formation of this wave regime manifests itself in, 
for example, fluctuations of the wave field at high frequencies and in the establishment of a 
near-Gaussian probability distribution function for the second sound wave amplitudes, as 
was observed in Ref . [44] . This regime is a close analogue of classical wave turbulence [46] . 
We refer to it as the second sound turbulence. 

The nonlinear dynamics of second sound involves the interactions of sound waves on 
different scales. Hence the full equations (7) should be solved numerically. In these simu- 
lations, we consider one-dimensional longitudinal second sound waves in a long cylindrical 
resonator [22]. We neglect the possible generation of transverse, radial modes. Note, how- 
ever, that the excitation of transverse modes in a cylindrical resonator can be important 
for non-planar, three dimensional waves [70]. The longitudinal resonant wave vector is 
k n = irn/L where L is the length of the resonator and n is the number of the resonant wave 
mode. The equations of motion (7) governing non-local energy balance in the second sound 
system read 



.db n = 



i — ^ t ^n,m,na (Pni bni &n— n\— raa 26ni6 n2 <5 ni _ n2 — n) "i^lnbn ~\~ Fd- (^) 



Here, V n ,m,n2 ot2{ nn i n 2) l l 2 describes the three-wave interaction, 7 n = vn 2 models the 
viscous damping of second sound, oc W is the amplitude of the force driving the n-th 
resonant mode [22, 41]. In the present model, the wave damping was taken into account at all 
frequencies, a feature that is of key importance [22] for a correct description of the formation 
of a cascade of nonlinear sound waves with increasing amplitude of the driving force. The 
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Figure 1: (a) Steady-state second sound power spectral amplitudes <5T W calculated numer- 
ically from Eq. (17) for four different driving force amplitudes: Fd = 0.01 (triangles); 0.05 
(diamonds); 0.1 (circles); and 0.3 (squares). The dashed line corresponds to ST^ oc uj^ 1 . 
The arrow marks the boundary frequency Wfe/2-7r of the inertial range for Fj = 0.3. The 
inset shows the calculated dependence on Fj of the standing wave amplitude 5T: solid line 
- nonlinear waves, 02 < (T > T a ); dashed line - linear waves, ai = (T = T a ). (From 
Ref. [22].) (b) Formation of the inverse cascade of second sound waves spectrum computed 
for t = 1000. Main plot: the wave spectrum in presence of the inverse cascade computed 
for N = 10 modes in the resonator. The vertical arrow indicates the frequency u?(\rive 

of the 

driving force. Inset: oscillations of temperature at a wall of the resonator. The horizontal 
arrow shows the driving force period. 

wave spectrum is calculated as 5T U = B n (b n + b* n ) where B n = (oj n / L^idS/dT) 1 / 2 . We 
use the dispersion relation for linear waves, u n = U2k n , to compute the resonant frequencies 
of the cavity. 

Fig. 1(a) shows the evolution of the steady-state spectrum with increasing driving force 
amplitude F^, calculated for a.2 < (T > T a ). We used a periodic driving force of fre- 
quency equal to a resonant frequency of the resonator corresponding to the conditions of 
the measurements, taken as Wdrive/27r = 3000 in dimensionless units. The effective viscosity 
coefficient v was fitted to the measured value of the quality factor Q ~ 3 x 10 3 of the 
resonator, to facilitate comparison of the model results with those from the experiments 
[22]. Points on the plot correspond to the amplitudes of the peaks in the spectrum. It is 
seen that, at small driving amplitude ~ 0.01 (triangles), viscous damping prevails at all 
frequencies and a turbulent cascade is not formed: the amplitude of the second harmonic 
is an order less that the amplitude of the main harmonic. In this regime the wave shape 
is close to linear. At intermediate driving amplitude ~ 0.05 (diamonds) nonlinearity 
starts to play a role at frequencies of the order of driving frequency, and a few harmonics 
are generated. At high driving amplitudes F^ > 0.1 (circles and squares) a well-developed 
cascade of second sound waves is formed up to frequencies 30 times higher than the driving 
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Figure 2: (a) Evolution with time t of the spectral amplitude during the build-up of second 
sound turbulence, (b) Points: dependencies on t of the spectral amplitude at driving 
frequency, 5T\, and of higher harmonics, 5T n with n > 1. Lines: dependencies 5T n oc t 2n . 
The computations were for i 7 ^ = 1. From Ref. [45] 



frequency. 

2.3.2 The inverse cascade 

To capture the main characteristic features of the development of the inverse cascade while, 
at the same time, keeping the system simple enough to be analyzed in detail, we carry out 
a simulation within the framework of equations (17) taking account of the relatively small 
numbers of waves N: the summation on the r.h.s. of Eq. (17) was made for m, n<i ranging 
from 1 to N. In what follows we present the results obtained for N = 10. We assume that 
a periodic driving force is applied at the 6th resonant frequency. 

It is was observed in simulations that, for sufficiently high driving amplitude W, low 
frequency harmonics at uj < ^drivc are formed in the wave spectrum (see Fig. 1 (b)). Con- 
sistent with the discussion above, this is a manifestation of the formation of an inverse 
cascade in a superfluid. The numerical results describe very well the data obtained in the 
experiments on the inverse cascade of second sound waves in He-II (see Sec. 2.4.3 below). 

It is clearly evident from Fig. 1(b) that the inverse cascade is responsible for the forma- 
tion of high-amplitude, low-frequency, subharmonics (cf. the experimental results shown in 
Fig. 4(b)). As it is shown in Refs. [23, 24], the inverse cascade develops through formation of 
isolated low-frequency waves of higher amplitude than is typical of the waves around them. 
These higher-amplitude lone waves can be considered as the acoustic analogue of the giant 
"rogue" waves that occasionally appear on the ocean and endanger shipping. Their origin 
lies in a decay instability of the periodic wave, i.e. a similar mechanism to that proposed 
[71, 72] (modulation instability) to account for the creation of oceanic rogue waves [73]. 
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2.3.3 Turbulence formation 



To capture the dependence on time of the formation of the direct cascade, we compute 
the evolution of the second sound wave spectrum with time for the driving force F = 
1 (in numerical units). The steady-state spectrum established at large times, i.e., after 
all transient processes are finished at t > tf orm ~ 5, is shown in Fig. 2(a) by circles. 
The numerical unit of time we use in simulations is equal, in order of magnitude, to the 
characteristic time T2 taken at the driving frequency. Figure 2(a) also demonstrates the 
evolution of the spectrum during the build-up process. Figure 2(b) shows the dependence 
on time of the amplitude of the wave, 5T n at the driving frequency (n = 1) and of two higher 
harmonics at n = 2 and 3, computed for the evolution presented in Fig. 2(a). (We recall that 
the second sound oscillations observed in experiments are fluctuations in the temperature 
of the superfluid). It is clearly evident that the dependence of the wave amplitudes on time 
during the formation process t < 2 is well described by a power-law function 5T n oc t m , 
where m = 2n — 1. Note that this power-law dependence has also been found in analytical 
computations through a small-time expansion in Eqs. 17 [45]. 

The nonstationary process through which second sound turbulence is formed has recently 
been investigated in experiments Ref. [24]. It was found that the turbulence build-up 
process is well described by the self-similar dependence (15) and is in agreement with Fig. 
2. From a comparison of experimental data with Eq. (15) the exponent was estimated to 
be p = 5. 

2.4 Experimental observation of the direct and the inverse cascades 

2.4.1 Experimental arrangements 

Detailed experimental investigations of second sound turbulence in superfluid helium in a 
high-quality resonator have recently been undertaken [22,23,45,53-56]. In these studies, 
a thin-film heater was utilized to generate the second sound, and a thin-film bolometer (a 
fast-acting superconducting thermometer) was used as a detector: see Fig. 3. The use of 
a high-Q resonator enables one to create nonlinear standing second sound waves of high 
amplitude with only small heat production at the source. The resonator was formed by a 
cylindrical quartz tube of length L = 7 cm and inner diameter D = 1.5 cm. The film heater 
and bolometer were deposited on the surfaces of flat glass plates capping the ends of the 
tube. The heater was driven by an external sinusoidal voltage generator in the frequency 
range between 0.1 and 100 kHz: The frequency of the second sound (at twice the frequency 
of the voltage generator) was set equal to one of the longitudinal resonant frequencies. 

2.4.2 The direct cascade of second sound turbulence 

For a small ac heat flux density W < 1 mW/cm 2 , corresponding to a standing wave am- 
plitude less than 0.5 mK, a nearly linear regime of wave generation was observed: the 
amplitude 5T of the standing wave was proportional to the heat flux density W . Increase 
of the excitation above a few mW/cm 2 led to large deviations from the linear dependence 
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Figure 3: Schematic diagram illustrating the experimental arrangements. The cylindrical 
walls of the resonator are not shown. From Ref. [65] 



5T oc W, however, and to visible deformation of the initially sinusoidal standing wave, 
accounted for by the formation of multiple harmonics in its spectrum. Fig. 4(a) shows 
some typical experimental results at 2.079 K. The formation of multiple harmonics in the 
spectrum corresponds to the establishment of second sound turbulence in superfluid helium 
in the resonator. 

It is clearly evident from Fig. 4(a) that the main spectral peak lies at the driving 
frequency Wdrivo and that high-frequency peaks appear at its harmonics, uj n = thrive x n 
with n = 2,3,... It can be seen that a cascade of waves is formed over the frequency 
range up to 80 kHz, i.e. up to a frequency 25 times higher than the driving frequency. The 
dependence of peak height on frequency may be described by a power law <5T W = const x uj~ s 
at frequencies lower than some cutoff frequency ujb/2ir ~ 5.5 x 10 4 Hz. For sufficiently high 
heat flux densities W > 10 mW/cm 2 the scaling index tends to s ~ 1.5 ± 0.3. 

2.4.3 Build-up of the inverse cascade 

Formation of the inverse cascade of second sound turbulence was observed when a high 
enough periodic heat flux density was applied at a small detuning from the resonance. It 
was found that in this case turbulence formed in two stages. First, the direct cascade 
developed after the heat flux was switched on, as shown by the blue spectrum in Fig. 4(b). 
The green arrow on the main plot marks the driving frequency Wdrive- The moment t = 0.397 
s, at which the the driving voltage was switched on, is labeled by a black arrow in the inset. 
At the second stage, subharmonics of the driving frequency at oj < arrive were generated, as 
demonstrated by the orange spectrum on Fig. 4(b). Orange triangles in the inset show the 
dependence on time of the wave energy accumulated in the low- frequency region cj < Wdrive- 
The reduction of wave amplitude seen in the high frequency spectral domain (blue triangles) 
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is indicative of the onset of energy backflow towards lower frequencies, i.e. a sharing of 
the input energy between the direct and inverse cascades. The decrease in energy at high 
frequencies in 0.397s < t < 1.3 s is attributable to relaxation processes in the direct cascade. 
Redistribution of wave energy due to sharing of the energy flux between the direct and 
inverse cascades starts at t = 1.3 s. The results of our numerical computations discussed 
above in Sec. 2.3 are in agreement with these experimental observations. 




Figure 4: (a) Power spectral amplitudes 5T W of standing waves recorded at T = 2.079 K 
when driving at the 31st resonant frequency, Wdrive/2vr = 3130 Hz. The ac heat flux density 
from the heater was 22 mW/cm 2 . The dashed line corresponds to <5T W oc u; -1 ' 5 . The arrows 
indicate positions of the viscous cutoff frequency, (b) Formation of the inverse cascade 
of second sound turbulence. The lower blue spectrum shows the direct cascade only; the 
upper orange spectrum shows both the direct and inverse cascades formed at later stage of 
the turbulence formation. The green arrow indicates the fundamental peak at the driving 
frequency. Inset: evolution of the wave energy in the low- frequency and high frequency 
domains is shown by the orange squares and blue triangles respectively; black arrow marks 
the moment of time at which the external drive was switched on. From Refs. [22, 23] 



3 Bose-Einstein condensation and superfluidity of polaritons 
in a microcavity 

3.1 Microcavity polaritons 

In this Section, we consider another important example of a quantum system, in which the 
Bose-Einstein condensation occurs at low enough temperatures - the exciton polaritons in 
a quantum well (QW) embedded into an optical microcavity. The microcavity just consists 
of two mirrors, which are built opposite to each other in a semiconductor heterostructure 
[74-76], as shown in Fig. 5(a). The two mirrors are typically fabricated as sets of dielectric 
layers forming a Bragg reflector [74]. The cavity is open in the transverse direction, so 
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Figure 5: (a) Exciton in an optical microcavity: a quantum well is placed between two 
mirrors composed of multilayer Bragg reflectors, (b) Parabolic red dash curve: dispersion 
of photons in a planar microcavity, k\\ is a component of the photon wave vector parallel to 
the mirrors. Horizontal red dash line: Energy of excitons in a quantum well resonant with 
the cavity photon energy. Compared to the dispersion of the cavity photons, the exciton 
energy is essentially constant near k\\ = 0. Solid blue lines: upper and lower polariton 
modes formed by anticrossing of the photon and exciton states. 

that light is free to move out of the cavity if it propagates parallel to the mirrors, but it 
is confined to specific resonant modes in the longitudinal direction perpendicular to the 
mirrors. The quantum well is placed in between these two mirrors, as seen in Fig. 5(a). 
The excitons, which are the bound states of an electron in the conduction band coupled to 
a hole in the valence band, are created by external laser pumping [75]. Let us note that 
by designing Bragg reflectors of proper thickness, it is possible to tune the fundamental 
frequency of the cavity to be equal to the energy of the excitons; under these conditions 
light and excitons interact resonantly because of the non-zero electric dipole moment of 
the excitons, thus forming two new branches of excitations. These new branches, known 
as the upper and lower polaritons, are superpositions of the excitons and photons in the 
microcavity (see Fig. 5(b)) [75,77,78]. 

The effective mass of the lower polariton branch is given by the curvature of the band 
at zero wave vector in the longitudinal direction of the QW, i.e. at ku = 0, where k\\ 
is a component of the photon wave vector k parallel to the Bragg mirrors. Because of 
the small curvature of the band, the polaritons are extremely light relative to atoms and 
excitons (cf. discussion and references below). Under typical experimental conditions, the 
effective mass of a polariton can be as low as 10~ 4 x too, where too is the free electron 
mass [75]. Consequently, as we demonstrate below, the onset of BEC and superfluidity 
occurs in the system of polaritons at much higher temperatures than in systems of atoms 
or excitons [76, 79, 80]. 

It is worth noting that the quantum well in an optical microcavity behaves like a laser 
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because the light emitted from the polariton condensate is nearly coherent [76]. However, 
this special kind of laser does not require an inverse population of excitons and thus, the- 
oretically, it has no threshold [82]. Several recent experiments [14,15] have demonstrated 
spontaneous coherence in exciton-polariton systems in various two-dimensional (2D) semi- 
conductor microcavity structures. In addition, the resonant laser pumping can promote, in 
its turn, condensate formation as a result of nonlinear light scattering [81]. 

3.2 Bose-Einstein condensation and superfluidity of two-dimensional po- 
laritons in an in-plane harmonic potential 

The properties of polaritons have been studied in several theoretical works. The theory 
of polariton dynamics due to the polariton-polariton interaction has been developed in 
Refs. [83-85]. The crossover between lasing and polariton coherence has been studied in 
Refs. [86, 87]. Polaritons superfluidity and spontaneous linear polarization of the light emis- 
sion were predicted in Refs. [18, 88]. In these earlier studies, the coherent polaritonic phases 
were analyzed in unrestricted 2D system. 

At finite temperatures, T > 0, there is no BEC in an infinite two-dimensional (2D) 
system [59]. Nonetheless, a BEC phase transition can be obtained in a finite-sized system 
in the presence of a confining potential [89,90]. Recently, polariton systems in a harmonic 
potential trap have been studied experimentally in a GaAs/AlAs quantum well embedded 
in a GaAs/AlGaAs microcavity [91]. For this kind of trap, the exciton energy was tuned 
by applying stress to the host material, thus changing the inter-band gap energy. A BEC 
of polaritons in a quantum well has been observed experimentally in such systems [92]. 

Below, we briefly review the theory of Bose-Einstein condensation in the polariton sys- 
tem [16]. In the computations, we take advantage of the Hamiltonian approach. As it was 
pointed above, in a quantum system the Hamiltonian normal variables should be changed 
to operators that act in the occupation number space [59,93]. The operators of creation 
and annihilation of polariton in a lower energy band, ffifo) and V>(r), are defined as 

^(r) = 5> P e^, V f (r) = E^P e "^- ( 18 ) 
p P 

Here, we use an expansion over running plane waves, pp and p F that are the annihilation and 
creation operators of polaritons with the momentum P = hk. The latter are the quantum 
analogues of the normal coordinates and ajjl at the wave vector k. The operators V^(r) 
and ^>(r) taken at the same moment of time t satisfy the Bose commutation relations [93, 94]. 

As we note above, we consider a "cloud" of excitons trapped in a semiconductor potential 
well, described by the function V{r) with r being the distance from the centre of the 
potential. The external potential V(r) can arise, for example, from a shift of the exciton 
energy by application of inhomogeneous stress to the host material [91]. In our model, 
we also suppose that the photon states in the cavity are unaffected by the stress. In the 
computations, we use the simplest, harmonic approximation for the potential energy in the 
form V(r) = (l/2)7r 2 ; the parameter 7 characterizes the strength of the interaction. Below 
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we also refer to 7 as a "spring constant" , in analogy with the classical relation between the 
energy of a linear spring and the displacement r. 

In the case where that both the polariton momentum P and the spring constant 7 are 
small the effective Hamiltonian for the lower-band polaritons in the parabolic trap V(r) 
acquires the following form [16]: 

r - / n 2 A \ - r/ 0) 

H cS = J dr V f (r) ( + V eS (r)\ V(r) + -f- J dr ^ {v)^ {v)^{v)^{v), (19) 

where the effective mass of a polariton is defined as 

M-s 1 = \ (M- 1 + (c/n)\ c /(2nhj) . (20) 

Here c is the speed of light, n = y/e is the refraction index in the microcavity, e is the 
dielectric constant, and Ac is the length of the microcavity (i.e. the spacing between the 
Bragg mirrors), M = m e + mh is the exciton mass, m e and m/j are the effective masses of 
electron and hole. 

The effective external potential is equal to V e g(r) = ^jegr 2 with 7 e g- = 7/2. The last 
term in Eq. (19) describes the interaction of polaritons (a quantum analogue of the four- wave 
scattering considered in Sec. 2). The polariton-polariton interaction is purely repulsive, with 
the hard-core contact potential U c g(r — r') = U^S(r — r'), where U^g = 3e 2 a2£>/(2e) is 
the effective interaction constant, e is the charge of electron, a,2D = ^ 2 /(2m e _ft,e 2 ) is the 
two-dimensional Bohr radius of the exciton, and m e _h = m e mh/(m e + mh) is the exciton 
reduced mass. 

As can be seen from Eq. (19), the effective Hamiltonian for the trapped polaritons maps 
onto the Hamiltonian of a weakly-interacting 2D Bose gas in the confining parabolic trap 
V e ff(r). As it shown in Ref. [89], in a parabolic trap the Bose particles can form a Bose- 
Einstein condensate, and hence enter a superfluid state at temperatures below a critical 
temperature Tc°^ = fik^ \/2)^ c g N/itM e g , where ks is Boltzmann's constant. According to 
Eq. (19), this expression gives an estimate for the temperature of the BEC transition, 
for weakly interacting polaritons. 

In this approach, we can estimate the number of polaritons in the condensate fraction at 
T < T c (0) assuming that the polariton system is in thermal equilibrium. The approximation 
is valid if the characteristic time for polariton-polariton scattering is short compared to the 
polariton lifetime, in which case the polariton gas rapidly thermalizes and can therefore be 
characterized by well-defined thermodynamic temperature. Note that, this polariton tem- 
perature can be substantially higher than the lattice temperature of the host semiconductor 
sample. Although the polariton lifetime is typically short, ~ 2 ps, the polariton-polariton 
scattering time is even shorter at sufficiently high polariton density. A recent study [14] 
has demonstrated that, in an actual experimental system, this requirement is satisfied and 
polaritons are indeed characterized by a well-defined temperature. 

Below we restrict our consideration to the range of temperatures below but ef- 
fectively above the lowest energy level in the confining potential, T S> huo/kB, where 



19 




Figure 6: (a) Condensate fraction Nq/N of polaritons in a parabolic trap in an optical 
microcavity as a function of temperature T, calculated for three values of trapping potential 
strength, (solid line) 7 = 760 eV/cm 2 ; (dotted line) 860 eV/cm 2 and (dashed line) 960 
eV/cm 2 . (b) Ratio of BEC critical temperature to the square root of the total number of 
magnetopolaritons T^ 1 / \/~N as a function of magnetic field B at different trapping potential 
strength 7. We assume that the environment around the graphene layer is GaAs with 
e = 12.9. From Refs. [16, 17] 

w o — (7cff/-^eff) 1//2 - For example, by taking 7 = 100 eV/cm 2 as an estimate for typical 
experimental conditions, this implies T S> 0.3 K. The latter condition allows us to neglect 
geometrical quantization of the polariton energies in the confining potential, and hence use 
the quasiclassical approach to the problem [16]. Figure 6(b) shows the fraction Nq(T)/N 
of polaritons in the condensate at a given temperature T for the experimental parameters 
provided in Ref. [16]. The computations were done for a weakly interacting Bose gas of 
polaritons at n(Q)a^ D <C 1, where n(0) is the density of polaritons at the center of the trap. 
In this figure it is clearly evident that a Bose-Einstein condensate exists in the polariton 
system up to T > 20 K, i.e. at temperatures much higher than the temperature of the 
normal-to-superfluid transition in 4 He, T\ = 2.17 K. 

3.3 Trapped polaritons in graphene in a microcavity in high magnetic 
fields 

Recent advances in the fabrication of carbon-based materials have enabled the production of 
graphene, that is, a two-dimensional honeycomb lattice of carbon atoms forming the basic 
planar structure [95-97]. Graphene is a semi-metal that is characterized by an unusual 
electronic structure. In particular, the electronic excitations with energy close to the Fermi 
energy (the last occupied energy level) have a linear dispersion with virtually zero effective 
mass. Therefore, these excitations in graphene represent a condensed-matter analogue of 
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Dirac fermions (like neutrinos) in standard quantum field theories [98]. Because of this, and 
due to high potential for the use in electronic devices, graphene has stimulated a considerable 
amount of theoretical and experimental study [98-100]. In this Section, we briefly describe 
the properties of trapped polaritons is a single graphene layer (GL) in an optical microcavity 
in a high magnetic field (i.e. magnetopolaritons: see below). In particular, we show that 
this system can undergo a superfluid transition at finite temperature [17]. 

It is known [101] that charged particles can readily be trapped in a potential well in an 
external magnetic field directed perpendicularly to the layer. Because of the quantization of 
particle motion in a magnetic field, the energy of the charged particles can only take discrete 
values, called Landau levels, at E n = (n + l/2)hcoB — 2fiBS z B, where B is the magnetic field 
strength, n is the level number, and ojb = \e\B/mc is the cyclotron frequency [101], m is the 
mass of the particle, \ib is the Bohr magneton, and s z is a projection of the electronic spin 
on the direction of the magnetic field. In a magnetic field, multiple electrons can occupy 
the same Landau level. The degeneracy of an individual Landau level i.e. the maximum 
number of electrons on it, is proportional to the magnetic field strength, v = \e\BS/2ir%c 
with S being the layer area [59, 101]. In full analogy with consideration given in Sec. 3.1 
above, the electrons and the holes in a graphene layer in magnetic field are coupled into 
excitons, which are referred to as magnetoexcitons [17]. 

In undoped graphene, in a perpendicular magnetic field (without an external electric 
field) and in the ground state, half of the 0th Landau level, i.e. v/2 states, is occupied by 
electrons, whereas all Landau levels at n > are empty. We label the Landau levels for 
electrons by numbers n > 0, and the Landau levels for the holes that are the positive charge 
carriers, by n < 0. Thus, the chemical potential of the system is positioned at the level 
n = 0. 

In the presence of an electric field, which is produced in experiments by applying the 
"gate" voltage to the system, the chemical potential of the magnetoexcitons can be changed 
in either of two ways: (a) an increase in the chemical potential in such a way that the 
resulting chemical potential is positioned between the 0th and 1st Landau levels; or (b) a 
decrease in the chemical potential so that it is positioned between the first negative (—1) and 
0th Landau levels [102]. In case (a) the magnetoexcitons are formed from electrons in the 
1st Landau level and holes in the 0th Landau level, whereas in case (b) the magnetoexcitons 
are formed from electrons in the 0th Landau level and holes in the Landau level at n = — 1. 
(Note that, by applying an appropriate gate voltage we can also use any other neighboring 
Landau levels at n and n + 1.) In both cases, all Landau levels below the chemical potential 
are completely occupied and all Landau levels above the chemical potential are completely 
empty. 

For a relatively high dielectric constant e of the microcavity, the magnetoexciton energy 
can be calculated by applying perturbation theory with the small parameter being the 
strength of the Coulomb electron-hole attraction [17]. The computation is similar to that 
[103] made for electrons and holes in a 2D quantum well in high magnetic fields. This 
approach is valid if the Coulomb electron-hole interaction energy in the single graphene 
layer, e 2 /(ers) is smaller than the energy difference between the Landau levels in graphene, 
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hvp/ri,, where r& = (hc/^B) 1 / 2 denotes the magnetic length of the magnetoexcitons in the 
magnetic field B, vf is the Fermi velocity, that is, the velocity of electrons of energy equal 
to the Fermi energy. This requirement is satisfied, for example, for graphene embedded in 
a GaAs microcavity, for which e ~ 13. (Note, that we restrict ourselves to consideration of 
high magnetic fields.) 

The effective Hamiltonian of the magnetopolaritons trapped by the potential V(r) in 
the graphene layer embedded in a microcavity is 



^eff = E — T^TT- + «V(T) , (21) 




where ft F and ftp are the lower polariton creation and annihilation operators (cf. Sec. 3.1). 
We see that the magnetic field strength only enters into the Hamiltonian (21) through the 
effective mass of the magnetopolaritons, which is 



^ (B) = 2 l < + TJ ■ (22) 

Here, Ac is the length of the cavity, and tub = 2 7 / 2 eh 2 /ir 1 / 2 e 2 rB is the magnetoexciton mass 
in graphene [17]. We suppose that the magnetic field is chosen so that the magnetoexciton 
and photon modes are in mutual resonance. Therefore, the length decreases with increasing 
magnetic field as Ac* oc B~ 1 / 2 . 

From Eq. (20) it follows that the effective magnetopolariton mass increases with 

the magnetic field as B 1 / 2 . The Hamiltonian (21) is similar to that for the pure polariton 
system (19). Meanwhile, it can be shown that in a high magnetic field the interaction 
between magnetopolaritons is negligibly small [17] and, because of that, that the fourth 
order term on the r.h.s. of Eq. (21) can be omitted. Thus the magnetopolaritons in a 
quantum well in a graphene layer form a virtually ideal quantum gas. 

As explained in Sec. 3.1, two-dimensional Bose particles in a trapping potential un- 
dergo Bose-Einstein condensation at some critical temperature T = T c . Computations 
similar to those made in Sec. 3.1 give the following expression for the condensate fraction 
of mag netopolaritons at T < T c (0) [17], 

N (T, B) = N- -± — 2 >- (k B T) 2 . (23) 

3n 7 

In Eq. (23) N is the total number of magnetopolaritons, and gs and glfi are the 
spin and momentum degeneracy factors for electrons and holes in graphene, respectively. 
At T = Tc°^ the condensate fraction vanishes. Setting iVo = in Eq. (23), we estimate the 
critical temperature for the ideal gas of magnetopolaritons as 



T&>(B) = ±( ^ 
c { } k B \8nM eS (B 



2 Ar \V2 

(24) 
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According to Eq. (24), the critical temperature T c decreases with magnetic field as B 1 / 4 
and increases with spring constant as 7 1/2 . Fig. 6 (right) demonstrates the dependence 
of the normalized critical temperature for the magnetopolariton condensate on external 
magnetic field. In the computations, we used gf^ = g^ 1 = g£^ = gi h ^ = 2. 

It is important to note that taking account of higher order terms in the Hamiltonian 
(21) (which are of the relative order of 1/e) results in magnetopolariton-magnetopolariton 
interactions and hence superfluidity in the magnetopolariton system. 

3.4 Drag effects in the system of electrons and microcavity polaritons 

Because the exciton consists of a coupled electron and hole, and hence has zero net charge, 
it can not be moved directly by application of an external electric field. However, an 
electrical current in a neighboring, spatially separated, quantum well can entrain excitons 
and therefore induce their motion [105]. In this subsection, we consider such drag effects 
for polaritons in a quantum well in a microcavity and propose possible experiment on the 
generation of collective polariton flow. Moreover, because the excitons are entangled with 
the cavity photons, polariton motion causes changes in the distribution of light emitted from 
the microcavity. In effect, we can direct light by application of voltage to the quantum well 
[104]. 

In this setup, we consider two neighboring quantum wells embedded in an optical micro- 
cavity as shown in Fig. 7(a). The "upper" quantum well (QW1) is occupied by 2D electron 
gas, and the "lower" quantum well (QW2) is occupied by the excitons created by laser 
pumping. It is important to note that the quantum wells are not identical: QW1 are doped 
so that there are excessive free negative charges, whereas in QW2 the number of negative 
and positive charges are equal to each other. We will consider the case of relatively low tem- 
perature where the occupation of the upper polariton branch is exponentially small, that is 
ksT <C Ti£Ir where Q# is Rabi splitting, which is equal to half the transition frequency from 
the upper to lower polariton branches at ku = (see Fig. 5(6)). At such temperatures the 
polaritons in QW2 form a superfluid [105] . In the computations, we omit small effects due 
to the nonresonant interaction of photons in the microcavity with electrons in the upper 
quantum well QW1. 

By applying an electric voltage, an electron current is induced in QW1 (its direction is 
shown in Fig. 7 by the arrow). This current drags the normal component of excitons in the 
neighboring QW2 due to polarization effects [104, 105]. 

This entrainment effect can be characterized by the drag coefficient of the polaritons, 
jp, introduced as follows [105] 

i p = -D p Vn p + j p E, (25) 

where i p is the mass current of the polaritons, D p the polariton diffusion coefficient in the 
QW2 plane, n p is the two-dimensional density of thermally activated quasiparticles, which 
form the normal component of the polariton superfluid, and E is the lateral electric field in 
the upper quantum well QW1. Note that the drag effect is fully reversible: the polaritonic 
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Figure 7: (a) The schematic diagram for the drag effects in quantum wells embedded in an 
optical microcavity. At low temperatures the polaritons in the lower quantum well QW2 are 
dragged by the electron current in the upper quantum well QW1. (b) The drag coefficient 
7 P in V~ 1 s~ 1 in the system of superfluid microcavity polaritons and electrons as a function 
of temperature T in K and total polariton density ritot in ?n~ 2 . The interwell separation is 
D = 20 nm. We used the parameters for the GaAs/GaAsAl quantum wells: m e = 0.07mo, 
m h = 0.15m , M = 0.24m , and e = 13. From Ref. [105] 

current formed in the quantum well QW2 will generate a non-zero current of electrons in 
QW1 [105]. 

The dependence of the drag coefficient -y p on the total density n to t of the temperature- 
activated quasiparticles in the polariton superfluid, and on the temperature T, is demon- 
strated in Fig. 7(b) [105]. It follows from the figure that the drag coefficient grows with 
increasing temperature T and slowly decreases with increasing density ntot- It can De 
demonstrated by detailed computations [105] that the dependence of the 7 P coefficient on 
the both parameters is defined mainly by variations in the density of the normal component 
of the superfuid. In particular, 7 P tends to zero in the limit of very low temperatures, at 
which the normal component vanishes. We have to note that T here represents an average 
temperature of the system; the exact, local temperature varies in the direction of the electric 
field due to changes in the concentration of the quasiparticles (i.e. due to non-zero gradient 
of the normal density in the polariton superfluid.) 

Generation of the lateral polaritonic current due to the drag effect causes crucial changes 
in the photoluminiscence spectrum, that is, in the angular distribution of the photons 
leaving the microcavity. We recall that the mirrors in the setup shown in Fig. 5(a) are 
semi-transparent. Fig. 8 demonstrates the angular dependence of the intensity F(9) of 
photons emitted by the excitonic superfluid. We defined 9 as an angle between the k vector 
of the emitted photon and the unit vector normal to the Bragg mirrors. As illustrated in 
Fig. 8, a sharp peak in the photon intensity F{0) caused by photoluminiscence from the 
superfluid component is positioned at 6 = in the cases of both zero (Fig. 8(a)) and finite 
(Fig. 8(b)) polariton current i p . Nevertheless, the broad background distribution in F(9), 
which is generated by photoluminiscence of the normal polaritonic component, is shifted 
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Figure 8: Proposed ex- 
periment: the manifesta- 
tion of polar iton drag ef- 
fect through change of the 
angular distribution of the 
photons escaping the opti- 
cal microcavity. (a) The 
angular distribution of the 
photons escaping the op- 
tical microcavity without 
drag, (b) The angular dis- 
tribution of the photons es- 
caping the optical micro- 
cavity in the presence of the 
drag effect. F{9) is the 
light intensity emitted from 
the microcavity at an an- 
gle 8 with respect to the di- 
rection normal to the Bragg 
mirrors. From Ref. [105] 

for finite polariton current compared to the case of zero polariton current. In other words, 
we can control the angle at which the light is emitted from the microcavity by generating 
an electronic current in the quantum well QW1. 

This effect opens up an opportunity for designing thresholdless tunable micro-lasers, 
which may be of importance for applications in modern quantum electronics [104, 105]. 
Besides, the control of photons and/or excitons by the exciton-electron drag can be used 
for studies of the properties of the polariton and exciton system and, in particular, of 
superfluidity in the system. 

3.5 Discussion 

We presented above the current status of theoretical research of BEC and superfluidity 
of trapped quantum well polaritons in a microcavity. Besides, we report BEC of trapped 
magnetoexciton polaritons in GL and QW embedded in an optical microcavity in high 
magnetic field. In the both cases the polaritons are considered in a harmonic potential trap. 
The effective Hamiltonian of QW and GL polaritons in a microcavity in a high magnetic 
field, and the BEC temperature as functions of magnetic field, are obtained. It is shown 
that the effective magnetic mass of the magnetoexciton polariton depends on magnetic field. 

We have also discussed drag effects in the system of spatially separated electrons and 
excitons in coupled quantum wells embedded in an optical microcavity were discussed, and 
the corresponding experiments. 
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We considered how, for low temperatures, the Hamiltonian of the 2D exciton polaritons 
in a slowly varying external parabolic potential acts on the exciton energy and brings it 
into resonance with a cavity photon mode, corresponding directly to the case of a weakly 
interacting Bose gas with an effective mass and effective pairwise interaction in a harmonic 
potential trap. The condensate fraction is a decreasing function of temperature, as expected, 
and an increasing function of the curvature of the parabolic potential. Mixing with the 
photon states leads to a smaller lifetime for high-energy states, that is, an evaporative 
cooling effect, but it does not fundamentally prevent condensation. Since harmonic potential 
traps are now possible for microcavity polaritons, it should be possible to compare these 
calculations with experimental results for the critical density and spatial profile of the 
polariton condensate. 

Since microcavity polaritons form a BEC, light emitted from the condensate due to 
electron- hole pair recombination is nearly coherent. This coherence of the emitted light 
means that polariton BEC in a microcavity can be used in technology as polariton laser, 
since photon energy is pumped in incoherently, putting electrons into the excited states, 
and coherent light is emitted [75,76]. Moreover, because the emitted light frequency for 
the graphene-based polariton laser depends on the external magnetic field strength, magne- 
topolariton BEC in graphene in microcavity can be used for the design of a polariton laser 
tunable by magnetic field. 

4 Conclusions 

In conclusion, we have presented an introductory review of the non-equilibrium dynamics 
of quantum systems, a fast developing field. We devise an approach to the description 
of the out-of-equilibrium quantum liquid based on a Hamiltonian formalism of superfluid 
hydrodynamics. Through analysis of the quasiclassical kinetic equations for the correlation 
functions for the first and second sound wave modes we show that, in the non-equilibrium 
state i.e. where there is a constant energy input into the superfluid, the time evolution of 
the system is affected by nonlinearity. As a result of nonlinear wave interaction, the direct 
and inverse cascades corresponding to constant fluxes of the integrals of motion (the energy 
and the wave action), are established in the bulk superfluid. The stationary solutions of 
the kinetic equations that describe this state are similar to those that appear in classical 
nonlinear wave dynamics. However, in contrast to the classical picture, both the direct 
and inverse cascades appear as a result of three-wave interactions between sound modes. 
We demonstrate that analytic and numerical approaches describe well the results of ex- 
periments with nonlinear waves in superfluid 4 He. We also analyze the conditions under 
which Bose-Einstein condensate forms in the system of interacting exciton polaritons in 
optical microcavities in semicondictors, and in a graphene layer. In this work, the occupa- 
tion numbers of out-of-condensate quantum states are not large and the effect of quantum 
fluctuations is important; we therefore use an operator formulation of the Hamiltonian ap- 
proach. We show that because of small effective mass of the polaritons, the Bose-Einstein 
condensate in this system can form at a temperature much higher than the superfluid tran- 
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sition temperature in bulk liquid helium. We also studied the effect of an external magnetic 
field on the BEC condensation. The results obtained in this part are directly applicable to 
the systems and materials that are potentially useful for the design of novel microelectronic 
and of micro-lasing devices. 
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